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Abstract 

We study the azimuthal asymmetry of the unpolarized Drell-Yan in the framework 
of the T-odd functions. We find, on the basis of quite general arguments, that for 
|q_i_| «Q such an asymmetry decreases as Q~ 2 , where q_i_ and Q are respectively 
the transverse momentum and the center-of-mass energy of the muon pair. The 
experimental results support this conclusion. 
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1 Introduction 



As is well-known, unpolarized Drell-Yan (DY) presents an azimuthal asymmetry. 
This was seen, for example, in reactions of the type [1-3] 

7T~ N -> /i + /i"A, (1) 

where N is an unpolarized tungsten or deuterium target, which scatters a negative 
pion beam. Such an asymmetry was originally explained as a first order QCD correc- 
tion effect [1]. Recently, however, it has been attributed[4, 5] to the quark polarization 
in unpolarized (or spinless) hadrons[6, 7]. Such a polarization reads 

n = ^f/^, P i). (2) 

Here x and p_i_ are, respectively, the Bjorken variable and the transverse momen- 
tum of the quark, P is the hadron momentum and a Lorentz invariant function 
introduced by Mulders and Tangerman[8] (MT). Lastly /io is a factor with the di- 
mensions of a momentum, introduced for dimensional reasons and to be discussed 
below. Notice that this polarization, caused by a T-odd interference between two 
different amplitudes contributing to the process (1), is compatible with parity and 
time reversal invariance[6, 7, 9]. We limit our study to DY events such that 

l<ul«Q, (3) 

where q_i_ and Q are, respectively, the transverse momentum and the center-of-mass 
energy of the muon pair. A problem we like to focus is the Q 2 dependence of the 
above azimuthal asymmetry. This behavior depends crucially on the factor /i . MT 
assume it to be equal to the hadron rest mass, therefore they predict a substantial 
energy independence of the asymmetry. On the contrary, in this note we show that 
there are compelling reasons for defining differently that parameter; this implies that 
the azimuthal asymmetry decreases as Q~ 2 . As we shall see, data support such a 
conclusion. 

The paper is organized as follows. In sect. 2 we give the general formulae for 
the DY cross section, illustrating the parameters which show deviations from the 
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naive parton model. Sect. 3 is dedicated to the quark correlation matrix, containing 
"soft" information for the inclusive cross sections at high momentum transfers. In 
sect. 4 we determine the parameter /x , either by comparison between the correlation 
matrix and the density matrix, or starting from a simple model. In sect. 5 we 
calculate the azimuthal asymmetry parameter v (see formula (4) below) in terms of 
the polarizations II of the active quarks. Using the result of sect. 4, we obtain the Q 2 
dependence of this asymmetry, which we compare with experimental results. Lastly 
we draw a short conclusion. 



2 The unpolarized Drell-Yan cross section 

The DY angular differential cross section is conventionally expressed as 
X dcx 3 1 1 

— — = — (1 + \cos 2 9 + fisin29cos(j) H — usin 2 9cos2(j)). (4) 

a ail Air A + 3 2 

Here Q = (9, <fi), where 9 and <fi are respectively the polar and azimuthal angle of the 
momentum of the positive muon in a given center-of-mass frame of the muon pair. 
In particular, in the present paper we adopt the Collins-Soper (CS) frame [10], whose 
z-axis is along the bisector of the beam momentum and of the target momentum, 
while the x-axis is parallel to qj_. Moreover A, // and v are parameters, which are 
functions of the overall center-of-mass energy squared s, of q^_, of Q and of the 
longitudinal fractional momentum xf of the muon pair. On the left- hand- side, da/dQ 
is a shorthand notation for 

da da 

dn dQdQ 2 dx F d 2 q ± ' ^ ' 

where xp is the Feynman longitudinal fractional momentum of the muon pair with 
respect to the initial beam. In the naive DY model, where the parton transverse 
momentum and QCD corrections are neglected, one has A = 1, /j, — v — 0. Therefore 
deviations of such parameters from the above naive predictions can be attributed to 
a nonvanishing transverse momentum of the partons inside the hadron or to gluon 
effects. In particular A ^ 1 is due uniquely to the former cause. Since such deviations 
have been measured for all three parameters, it appears most natural to interpret 
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them, as far as possible, as an effect of the quark transverse momentum, essential in 
the polarization (2). In particular, such an interpretation of the azimuthal asymmetry 
accounts rather well [4] for the observed large size of v and small size of 2]. 



3 The quark correlation matrix 

The correlation matrix is a very important tool for calculating inclusive reaction cross 
sections at high momentum transfers. Indeeed, in the framework of factorization [11], 
this matrix contains all information concerning the "soft" functions entering the the- 
oretical expressions of such cross sections. For example, the DY cross section, which 
we are interested in, has the following convolutive expression: 

^ = \ — L^W v . (6) 

Here a is the fine structure constant, s the squared energy in the overall center-of- 
mass system and x a and x b the longitudinal fractional momenta of the active partons, 
with x F = x a — x b . Moreover L^ v and W^ u are respectively the leptonic and hadronic 
tensor, i. e., 

= + k v T - g^k ■ k, (7) 
k and k being the four-momenta of the two leptons, and 

W " U = IJ d2 P± Tr l®A(x a , P±h^ B (x b , q ± - P± h1 • (8) 

$^4 and <3>s are the correlation matrices of the active quark and antiquark, which 
belong to the two initial hadrons, denoted here respectively as A and B. * The corre- 
lation matrix may be parametrized according to the Dirac algebra, taking into account 
Lorentz and parity invariance. On the contrary, we do not consider restrictions due 
to time reversal, since we admit the so-called T-odd functions, i. e., some particular 
interference terms, which, without violating any symmetry, change their sign under 
that transformation. Obviously, not all the Dirac components will contribute to a 
single process. 

*For the sake of brevity, we omit the contribution of the annihilation between an antiquark of 
the hadron A and a quark of the hadron B. However this does not affect our conclusion. 
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In order to determine /xo, it is necessary to consider a polarized nucleon, e. g., a 
transversely polarized one. In this case the correlation matrix is parametrized as 



$ = $ ± ~ $ Qa _|_ $ Qfe + $ x + $ , 



(9) 



where 



Oa 



^P + [fli+ + ^±9lTl5fl+ + -/llT75[^,^H 



+ 



$06 
$1 



+ 



4v/2 

^ (ft + \±qH) h 



= -M{e + g Tl5 fi). 



1 

2/i 



M • c 



(10) 

(11) 
(12) 

(13) 



Here we have used the MT notations for the "soft" functions. Moreover we have set 

P 

A± = -S ■ p±/li , P±=P-P-PT75' ( 14 ) 



M being the nucleon rest mass, P and p the four-momenta of the nucleon and of 
the quark respectively, and S the Pauli-Lubanski (PL) four-vector of the nucleon. 
Taking the z-axis along the nucleon momentum, it results P = (M/2P+, P + , 0j_), p 
= (p~ ,p + , p_i_)) and pj_ = (0,0, px). n± are lightlike vectors, such that n + ■ n_ = 1 
and whose space components are directed along the nucleon momentum. 

The term <3> is T-odd. In particular, for an unpolarized nucleon, this reduces to 



(15) 



which corresponds to the quark polarization II (eq. (2)), as it is straightforward to 
check[6, 7]. It is important to stress that this polarization is due to a coherence effect, 
for example to one gluon exchange between the hadron A and the active parton of 
the hadron B, or vice-versa[12]. This effect can be factorized only if the condition (3) 
is fulfilled[4]. Moreover the factor /x must be the same for all the functions involved, 



in order to normalize them appropriately. This factor was assumed to be equal to 
the rest mass of the hadron[8, 13]. As we shall see in a moment, this is not the most 
suitable choice, in order to interpret the functions contained in the correlation matrix 
as probability densities or interference terms. 

4 Determining /i 

We follow two different procedures for deriving the correct expression of the parameter 
jiQ. The first procedure consists in comparing the correlation matrix with the density 
matrix, to which <£> reduces for non-interacting quarks. The second procedure is based 
on a simple, but rather general, model, without any ad hoc assumptions. 

4.1 The density matrix 

The density matrix of a free, on-shell quark in a transversely polarized nucleon reads 



its spin aligned with (T = 1/2) or opposite to (T = -1/2) the proton spin. For the 
sake of simpolicity, although not essential in our formulae, we consider a reference 
frame - called P-frame from now on - where the nucleon momentum has a modulus 
|P| = V >> M. In appendix we show that in such a frame one has 



P±= E <fr(z,p±)^ + m)(l + 2T 75 &). 



(16) 



T=±l/2 



Here m is the rest mass of the quark, such that p 2 = m 2 . 2TS q is the quark PL 
vector, with S 2 = —1. qx(x, p±) is the probability density of finding a quark with 




Here we have set 



9(z,pl)= Yl Qt(x,p±), 



5q ± (x,p ± )= ]T 2Tg T (x,p ± ), (18) 



T=±l/2 



T=±l/2 



C 1 = E q ±\fi' + ,fi'_] 8 in0', 



(19) 
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C2 = + ^ (l - - ft. + ^bW-]} «n£K (20) 

and 

sm0' = -p ± -S/|p|. (21) 

Lastly n' ± are defined analogously to n±, but with the space component along the 
quark momentum. 

4.2 Comparison with the correlation matrix 

We equate the coefficients of the independent Dirac operators in eqs. (17) and (9), 
taking into account the relation 

p = V2xVn + + P± + (V- 1 ) . (22) 

We get, for a free, on-shell quark[14], 

fi = fi=Q, \ ± h^ = sine'5q ± , (23) 

\±h^ T = (1 — €i)sin9'5q± : Ho^ihr = (1 — e-\)sin0' E q 5q±, (24) 

A±£?it = (1 - e 2 )sin6'5q±, X±g^ = (1 - e 3 )sin9'5q±. (25) 

Here e\ = m/E q , e 2 = m/xV and €3 = m/2|p| are the correction terms to the chiral 
limit, which are generally small for light quarks. The terms of order O [{m 2 + PjJ/P 2 ] 
have been neglected. As regards /x , we require the various functions to be normalized, 
in the chiral limit, like 5q±, which is a difference between two probability densities. 
Therefore we assume 

\ ± = sin9', (26) 
which, according to eqs. (21) and (14), implies 

A*o = |pj- (27) 

Two remarks are in order about iiq. 

• This parameter is frame dependent, as well as the correlation matrix and the 
density matrix; however, in a specific reaction, it can be expressed in terms of 
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Lorentz invariant quantities. In particular, for DY, in the center-of-mass system 
we have 

0o * f ■ (28) 

• /i is independent of the dynamics, therefore the presence of interactions among 
partons does not affect our result. 

4.3 A model for hj; 

A nucleon may be viewed as a bound state of the active quark with a set X of 
spectator partons. In order to take into account coherence effects, we project the 
bound state onto scattering states with a fixed third component of the total angular 
momentum with respect to the nucleon momentum, J z , and with a spin component 
s = ±1/2 of the quark along a given (axial) unit vector s. For the sake of simplicity, 
following Brodsky et al.[12, 15] (see also [16, 17, 18]), we assume X to have spin zero, 
moreover we choose a state with J z = 1/2. Then 

\J Z = 1/2; s; X) = a\ L z = 0; s; X) + (3\ <-, L z = 1; s; X). (29) 

Here — > («— ) and L z denote the components along P, respectively, of the quark spin 
and orbital angular momentum, while a and (3 are Clebsch-Gordan coefficients. Then 
the probability of finding a quark with J z = 1/2 and spin component s along s, in a 
longitudinally polarized proton with a positive helicity, is 

|(P,A=l/2|J 2 = l/2;s;X)| 2 = a 2 \(P, A = 1/2 1 — L z = 0; s; X)\ 2 

+ P 2 \(P,A=1/2\^,L z = 1;s;X)\ 2 + I, (30) 

/ = 2a[3Re [(P, A = 1/2 1 — >, L 2 = 0; s; *)<(<-, L 2 = 1; s; X) |P, A = 1/2)] . (31) 

We are especially interested in the interference term 7, which may be interpreted 
as a quark polarization, independent of the nucleon polarization. Indeed, eq. (30) 
implies that 21 is the difference between the probability of finding a quark with spin 
component s along s and the probability that the spin of a quark along the same 
quantization axis be —s. After partial wave expansion, I reads 

oo 

/ = 2 J2 Re 

i,v =o 
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ie- i<p0 AiBf, P t (cos9 )P^ (cos6 ). 



(32) 



Here Ai and Bi are related to partial wave amplitudes; moreover 9 Q and 0o are re- 
spectively the polar and the azimuthal angle of the quark momentum, assuming the 
polar axis along P and, as the azimuthal plane, the one through P and s. In the 
P-frame one has 

P^cosOo) ~ 1, P^cosOo) ~ (33) 

Then eq. (32) yields 

1 x (Ucos(f) + Vsiri(f) ) , (34) 



xV 

where U and V are real functions of x and pj_, made up with A[ and -B/. Since s is 
an axial vector, parity conservation implies U — 0. Therefore 

|p_ L | y 

' -Vsiri(p = =b— — |pj_ x s|, (35) 



xV xV 

where the ± sign depends on the sign of sm0o- Therefore the interference term I 
is T-odd. Moreover, setting s = P/|P|, and comparing eq. (35) with eq. (2), we 
identify with V and we get 

pL = xV. (36) 

Eqs. (36) predicts that the quark tranverse polarization in an unpolarized (or spinless) 
hadron decreases as V^ 1 . But in the center-of-mass frame, for Q » M one has 
xV ~ Q/2, therefore we recover the result (28). 



5 Azimuthal asymmetry 



Inserting formulae (7) to (13) into eq. (6), the differential cross section reads, under 
the condition (3), 

„2 



da 



1 



a 



dildQ 2 dxpd 2 q± (x a + Xb)s 12Q 2 ^ 



f 



H 



[1 + cos 2 6)F f + sin 2 6cos2(b 



(37) 



Here the sum runs over all the light flavors and antiflavors, u, d, s, u, d, s. Moreover, 
omitting the flavor indices, we have 



F = J d 2 p ± f 1A (x a , pD/ib^o, (q± - P±f], 



(38) 



where fiA and fis are the unpolarized densities of the active partons. Lastly 

H = J d 2 P± hi A (x a ,p 2 ± )hi B [x a , (<u - p ± ) 2 }S, (39) 

where 

S = [2p ± • n(<u - p ± ) • n - p ± • (q ± - p ± )] (40) 
and n = q_i_/|q_i_|. But eq. (39) implies that for qj_ — > H oc q^: 

H = iJ qi, (41) 

where H = H (x a , Xb, qi) assumes a finite value for qj_ = 0. Comparing eq. (37) 
with eq. (4) yields 

, = 254^. (42) 
Inserting eqs.(41) and (28) into eq. (42), we get 

v = A ^ = A p 2 , (43) 

Ao being a proportionality costant and p = |q_i_|/Q. Our prediction is compared with 
the experimental resuls of v at different energies, both as a function of Q at fixed 
|q_J << Q (figs. 1 and 2, where also the MT assumption is tested) and as a function 
of p (figs. 3 and 4). 

Comparison of our result with data confirms that the value of po deduced in sect. 
4 (see eqs. (27) and (28)) has to be preferred to the one assumed by MT. In particular, 
such a result allows a very simple interpretation of the behavior of v versus p, which 
would be quite difficult to fit with the MT ansatz. 
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Appendix 

Here we derive formula (17) for the density matrix of a quark. To this end we 
write formula (16): 

P±= E Qr(x,P±)l^ + m)(l + 2T 7 M, (A. 1) 

T=±l/2 

using the same notations as in subsect. 4.1. Moreover we define a quark rest frame, 
whose axes are parallel to those of the P-frame, defined in subsect. 4.1. Here the 
Pauli-Lubanski vector of the quark results to be S g = — S = (0, 0, 1, 0). Decom- 
posing S^°) = S into a transverse and a longitudinal component with respect to the 
quark momentum, we get 

= S = Z ± cos6' + usinO'. (A. 2) 

Here 

sind' = sin9 sirup, sind = -j-^-, sirup = f - — , (A. 3) 

IpI Ip±I 

P - (p±,xV), ^ = (0,t), E ± = (0,n), (A. 4) 

t = (sin9cos(p, sin9 sirup, cos9), (A. 5) 

n = (cos9cos<p,cos9 sirup, —sinO). (A. 6) 

In order to calculate S q in the P-frame, we perform a boost along the quark 
momentum. This boost leaves £j_ invariant and transforms v into p/m, where 

p=(\p\,E q t), = yWp 2 . (A. 7) 

As a result we get 



S q = S + 

having defined 



m 



sin9', (A. 8) 



m 



5 = —= — n 



1 



Now we substitute eq. (A. 8) into eq. (A. 1), taking into account the definitions (A. 
4) and (A. 9) of v and 5, and exploiting the relations — j)$ = l/2[fi,p r \ — p ■ S , p ■ S = 
p± ■ S and v = ^(n f + — n'_). As a result we get eq. (17). 
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Figure 1: The behaviour of the asymmetry parameter v vs Q, at constant |q_i_| << <5 
and ^/s = 19.1 GeV. Data are taken from ref. 1 and fitted with formula v = B /Q 2 , 
B = 2.34 GeV 2 . The dashed line corresponds to the MT ansatz. 
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Figure 2: Same as fig. 1, y/s — 23.2 GeF, -B = 2.52 GeV 2 . Data are taken from ref. 
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Figure 3: The behaviour of the asymmetry parameter v vs p = |q_i_|/<5, — 19.1 
GeV. Data are taken from ref. 2 and fitted with formula v = A p, A = 1.35 
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Figure 4: Same as fig. 3, = 23.2 GeV. A = 1.17 
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